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Abstract 

We solve the equation of motion of boundary string field theory allowing generic 
boundary operators quadratic in X, and explore string theory non-perturbative vacua with 
massive state condensation. Using numerical analysis, a large number of new solutions 
are found. Their energies turn out to distribute densely in the range between the D-brane 
tension and the energy of the tachyon vacuum. We discuss an interpretation of these 
solutions as perturbative closed string states. From the cosmological point of view, the 
distribution of the energies can be regarded as the so-called landscape of string theory, 
as we have a vast number of non-perturbative string theory solutions including one with 
small vacuum energy. 
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1 Introduction 



As a non-perturbative formulation of open bosonic string, boundary string field theory 
(BSFT) [II |2] was proposed as well as cubic string field theory (CSFT) [3]. In general, so- 
lutions of string field theories are quite important as they would provide non-perturbative 
vacua of string theory, to look at the true capability of string theory. 

Recently, the multiple D-brane solutions, which have greater energies than the trivial 
vacuum, were proposed jU [5] in CSFT. It would have a significance equivalent to the 
proof of the original Sen's conjecture [71 13 [10] , since the D-brane creation is thought 
of as a necessary ingredient for a complete non-perturbative formulation of string theory. 
To climb up the SFT potential hill instead of rolling down the hill to get to the tachyon 
vacuum, it is indispensable to treat the string massive modes. 

After the construction of the analytic solution for tachyon condensation [11], various 
analytic solutions in CSFT have been found [121 EI21 [II]- in recent times, analytic forms 
of lump solutions [15j HE] and multiple D-brane solutions were proposed. In BSFT, as 
well, an analytic solution for tachyon condensation and lump solutions have been found 

PETJOEIIIE]. 

To solve the equation of motion of CSFT, we encounter the infinite-dimensional equa- 
tion, which is hard to solve. In fact, there are some subtleties of proposed solutions 
[Ml El [251 [261 [271 EE]- On the other hand, there is a consistent truncation scheme which 
reduces BSFT to a standard field theory with a finite number of fields. The BSFT ac- 
tion was constructed also for boundary interactions quadratic in the worldsheet field X, 
corresponding to a subset of massive modes of open string [20] . 

The purpose of this paper is to solve the equation of motion of the BSFT action for 
the quadratic boundary operators. In contrast to CSFT, only the tachyon field plays 
a significant role in the BSFT exact solution for tachyon condensation and the lump 
solutions such that the analysis is rather simple. For this reason, it is natural to expect 
that one may obtain a new class of solutions by involving some more boundary operators, 
aiming at new string vacua and a construction of a multiple-D-brane solution. 

We adopt the BSFT action for quadratic boundary interactions with arbitrary number 
of derivatives on the worldsheet given in [20], and solve the equation of motion numerically 
to find homogeneous static solutions. The condensation of the massive fields is taken care 
of to their all orders. So the solutions are non-perturbative ones at the classical level of 
SFT, in the same sense as for the non-perturbative tachyon vacuum solutions of the BSFT. 
We discover a large number of new solutions of BSFT. Interestingly, those energies turn 
out to be smaller than the D-brane energy. Our analysis strongly suggests the existence 
of an infinite number of solutions. 
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We also find that an approximately uniform distribution of the energies of the so- 
lutions, which suggests a relation to closed string excitations at the tachyon vacuum. 
Furthermore, from a cosmological point of view, the distribution of infinitely many so- 
lutions is reminiscent of the so-called string landscape. It is intriguing that a solution 
with any small energy may be possible in BSFT, to reconcile the cosmological constant 
problem. We also find solutions with a part of Lorentz symmetry broken, which would 
serve as a realization of the old idea of spontaneous Lorentz symmetry breaking (and 
CPT breaking) through SFT [2D 1221 123]. 

This paper is organized as follows. In the next section we review the derivation of the 
BSFT action and derive the potential for the tachyon field and massive fields associated 
with generic quadratic boundary interactions on the worldsheet. From the potential, we 
obtain the equations of motion and solve them numerically in secfJJ We show plots of 
numerical results as well, to show the energy distribution of the solutions. In secJU we 
present a possible interpretation of the solutions as relevance to closed string states, and 
study properties of the solutions. Finally, sec]5] is devoted to discussions. 



2 Review: the BSFT action 

We give a short review of boundary string field theory (BSFT) based on [U12] El In addition, 
we summarize the derivation of the BSFT action for quadratic boundary interactions 
following [20]. 



2.1 Generic formulation of BSFT 

The dynamical variables of BSFT are boundary coupling constants Aj associated with the 
boundary operators O l of ghost number 1. The BSFT action is given by the solution to 
the equation 

= --T 25 dO d6'(0\6){Q Bl O(6')}) x . (1) 

% J *-* 

Here T25 is the tension of the D25-brane, O = Yli^iO 1 an d Qb is the BRST charge. 
(• • • )a denotes the correlation function in the two dimensional field theory on a unit disk, 
described by a bulk world-sheet action Sb u ik with boundary interaction terms: 

Sbuik+^TA, rwx0). (2) 



3 See [29] 130] for relevant formulas and derivations. For the supersymmctric formulation, see for 
example [THl EU I32J 133] • Recent proposals of the formulation includes [31] ■ 
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Here 9 is the angle parametrizing the boundary of disk. V* is the vertex operator associ- 
ated with the open string state IV 1 ): 

|V*) = 6_ 1 |O i ), (3) 

where \0 % ) is as well associated with the boundary operator O l . One can show that 
the action S satisfies the Batalin-Vilkovisky (BV) master equation and so has a gauge 
symmetry [T]. 

It was given in [2] to write down the action directly for general V constructed only 
from matter fields. For general choice of V, 

oo 

{Q,0(e)} = J2cd n cF n (e), (4) 

n=l 

where F n are some matter operators. Hence the ghost correlation functions appearing in 
(PQ) are of the form (c(9)cd n c{9')) . Due to the form of the ghost correlation function 

(c(0)c(0')c(0")> = 2 (sin(0 - 9') + sm(0' - 9") + sin(0" - 9)) , (5) 

(CQ) can be written in the form of 

— = -T 25 dd d6> (V\6) (A(9 f ) + cos(0 - 9')B{d') + sin(0 - e')C{9'))) x , (6) 

% %) */ 

where A, B,C are linear combinations of F n . The operator A has the expansion in term 
of a basis {V(0)} of matter operators, 

A(9) = J2^V\9). (7) 

j 

Then the action is given by 




where Z is the partition function of the world-sheet theory (J2]) and g is a constant. 

2.2 BSFT action with generic quadratic boundary interactions 

Next, following Li and Witten [20], we derive the BSFT action for the most general 
quadratic boundary operators. Note that the quadratic part gives a free CFT on the 
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worldsheet, so the truncation of string theory to the one with generic quadratic bound- 
ary interactions is a consistent truncation [18]. Solutions of the BSFT with the generic 
quadratic boundary interactions amounts to solutions of the full theoryQ 
The generic quadratic operators ar^l 



O = c{6)V{6) , V(9) = Y + 4~ 1 : X ^ 9) r d9 ' u "^ 9 ~ d ') X ^') ■■ ■ (9) 
Here :: stands for the normal ordering as 

: X^6)X U {6') := X^9)X V {&) + 2^ 1b |1 - e l ( e -^\ . (10) 
We use the standard closed string action 

S'bulk = J d 2 aVh (^—h ab ^d a X,d b X v + ^b ab V a c b ^j . (11) 

Here, {a, u yw [Q — 6')} is a set of boundary couplings and h a b and r)^ = diag(— , +,+,.. .) 
are the metrics on the world-sheet and the target space respectively. Without loss of 
generality, we can assume u^ u {9) = u u ^{—9). Following the above terminology, {A,;} = 
{a,«f } with itf = uT k = d6u^(6)e- ike , and@ 

V a (6) = d a V = ^, V$ v = d uT V = ^ : X,{9) J^d6>e^ e -^X U (6>) : . (12) 

It is notable that the V^ v can be formally expressed as a linear combination of quadratic 
local boundary operators X^d r X v , which are the vertex operators corresponding to a 
constant field strength and to a set of massive modes of open string for r > 0. Since the 
world-sheet action is quadratic in X, we can solve this theory. The boundary condition 
is deformed by the boundary interaction as 



r-2-rr 

o = d r x»(8)+ / de'u" v {e-e')x v {e'), (13) 

Jo 

where r is the radial coordinate of the unit disk. The Green's function satisfying this 

4 This was the important observation for the proof of the Sen's conjecture [B] by BSFT [T71 [T51 IT5] . 

5 It is not necessary to normal order V^ v since it doesn't have singularity for generic u^ v {9) regular at 
9 = 0. In fact, in [20], V^ v was not normal ordered. The normal ordered form, however, is useful to see 
how the action reduces to the one associated with V = ^ + ^] u^T]^ : X^(6) :, so we adopt it 
in this paper. 

6 Until just before (|27p . we treat for all k £ Z as independent valuables when we take derivatives. 
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boundary condition is 

G^(z,w) = (X l ,(z,z)X u (w,w))x 

= --^Tlfiu (l n \ z - M 2 + In 1 1 — zw\ 2 ) - a'A 0y 



fjjyV 



-a'J2 {A k ^(zw) k + A„ k ^(zw) k ) , (14) 



fc=i 

where 



Here (l/(|/c|r/ + itfc)) Mi , is the inverse matrix of (|A;|?7^ + w^T)' 



i 

Notice that the correlation function in (Tl4l) is evaluated with the bulk action (llip with 
the boundary terms given by (JH])- The partition function of the world-sheet theory is 
determined from the differential equation 

Q 1-2-K />27T 



\nZ = --I7 / ^ / Ve lfe ( e - e ')(: X^X^ff) :}x = \ A k>1/fM . (16) 



In the last equation, we have used (1141) and A kfiu = A_ kvii . By integrating this differential 
equation, we obtain the partition function 

00 

Z = Mdet(u )- 1/2 e- a Y[ e fc " ltr fo- u *> det(l + k^r, ■ u k )- 1 . (17) 

k=l 

Here (rj ■ u k )^ = T]npU> p k and M is the normalization constant determined by demanding 
that the partition function reduces to V26, the volume of the target space, for a = u k = 0. 
The factor det(wo) -1 / 2 i n the partition function originates from the integration over the 
zero modes 

J d 26 xe~^<" x ^ = (2W) 13 det(M )" 1/2 . (18) 

Thus in the uq — > limit, the factor det(wo) _1/ ' 2 should be replaced by (27ra') _13 V26. This 
fact determines Af = (2na') 13 . It is interesting to see how (TIT]) reduces to the partition 
function with the boundary interaction 

25 
At=0 
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which was given in [2]. The vertex operator Vi is given by (Q by choosing u^ v {9) = 
u^rj^ u 5{6) or choosing vff = v^rf for all k. The so-called Weierstrass' product formula: 



= -e^ x Y[ (e k ~ lx (l + k^x)' 1 ) , (20) 

X k=l 

where 7 is Euler's constant, leads to 

25 

Me- a \\ (Vtfe 7/ r(/)j . (21) 



Z\ u v (e)=u^rf v 5(d) 



At=0 



This is nothing but the partition function for (|T9|) given in [2]. 

The remaining task to derive the BSFT action is to find {a{\ = {a a , off} and g. By 
applying {Q B ,X^{9)} = cd e X^9), one gets 

1 r 2n 

{QB,O{0)} = cd e cV{6)- — :cX li (6) / d9'u^(9-9')cd e >X l/ (9'): . (22) 

J 

Substituting this into (CQ) and using the ghost correlation function (j5j), we obtain 

A{d) = -V{9) - — *— : XJ6) [ d6' sin(# - &)vT(0 - e')d e >XJ6') : . (23) 
Ana' J 

This implies 

a a = -a , off = -k [vff +1 - - vff . (24) 

g can be determined as follows. According to (151), 



Jh. r d g {A (9')) = T 25 V a, JL Z , (25) 

27r Jo J d\i 



QQ* p2tt p'2ty 

da 



where we have used (JTj) and the fact that ^-integrals of cos(# — 9') and sin(# — 9') vanish. 
On the other hand, using d a Z = — Z and (I24p . the derivative of OH]) with respect to a is 

Consequently, we obtain (7 = 1 and 



S = —T< 



25 



( a + 1 + E K+i - - «sr) z • ( 27 ) 
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Since we are interested in homogeneous static solutions, we remove — (uQ l> /2)dZ/du^ 1 ', 
which leads to the kinetic term of the tachyon field T, from (127|) . Hence the potential 
term is 



U = 
where 



// oo \ oo 
d 26 x e' T [T + 1 + /T^r II e fe " ltr(r? - Ufc) det(l + k^rj ■ u,)- 1 , (28) 
\ 1 1 / i, i 



k=l I k=l 



P? = \< - < , M>2 = < = \k « +1 - - < ■ (29) 

Setting T = a + u lu x fJi x l ,/2a' , one can reproduce (1271) except the kinetic term. In partic- 
ular, —UQ V dZ/duQ V in (127)1 is obtained from T in the parenthesis in the potential (|28|) . 

We further restrict our attention in the the case where = u^jf" . Since the non- 
diagonal parts of tiff always accompany the other non-diagonal elements of u^" , this 
restriction is consistent in the sense that dV/du^ u = for \x ^ v. The potential for the 
tachyon field T and the diagonal elements of Uk is 

(25 oo j . % \ oo 25 

r+i-EE « [uh? t ~ i) II II + ^'O- ■ (30) 
— k — 1 / k — 1 [J* — 

where 

X i 

&\ = 2 U 2-<^ Pk>2 = 2 k K+l - U k~l) - U k ■ ( 31 ) 

Here we focused on the homogeneous static fields and performed the integration f d 26 x. 



3 The solutions of the equations of motion 

In this section, we solve the equations of motion derived from the non-perturbative po- 
tential ( )30l) of the BSFT. Since we have infinitely many degrees of freedom, we adopt 
an approximation and solve them numerically. We find a large number of solutions, and 
those solutions have peculiar properties. First, we present the equations of motion, and 
then show how to solve them numerically with an estimate of the validity of the approx- 
imation. Then finally we study the peculiar properties of the energy distribution of the 
solutions. 



3.1 The equation of motion 

To find solutions of the equations of motion with non- vanishing Uk, we first solve dU/dT = 
0. This gives two solutions, the first one is 

25 oo , , 

^ee4^4)- < 32 > 



H=0 k=l 
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and the second one is T = oo with arbitrary^ Obviously the second one is the 
tachyon vacuum solution [18] since the potential energy ( 130 1 ) vanishes, and here we see 
the consistency with the truncated solution of [IB] explicitlyjj 

In the following we consider only the first solution (132|) . Substituting this solution 
back into (13"0"j) gives 



U = T 25 V W e - 



(33) 



where 



/(«) 



g( ft (rr^4)- rV+log ( 1+rl "*») 

5". + E(ft-(^)+ l0 6( 1 + *^ 

fc=l 7 



and 



ft 



-U2-«l, /?fc>2 = ^ («fc+l 



(34) 



(35) 



2 ' - 2 

Here, m in f(u) is the shorthand notation of {«*;}. In the second line of (1341) . we have 
used Y^h=i( u k+i ~ u k~i) = —u\. Clearly, solutions of the equations of motion derived 
from the potential U are the stationary points of X]^=o Now the potential energy 

at a stationary point is in the form of 



U* = T 25 V 26 e"^M=o/K ) 



(36) 



where is a complete set of solutions of df /du = labeled by n. 

Our next task is to solve df /du = 0. Since this is an infinite dimensional equation, it 
is difficult to solve it analytically. However, we first note that we find a solution 



u 







(37) 



for all k > 1 and This implies T = by using ( 152|) . so, the solution is nothing but 
the trivial vacuum of the original D25-brane. It is important that the trivial D25-brane 
solution and the tachyon vacuum solution are allowed in our generalized scheme, as a 
check of the consistent truncation of the BSFT. 



7 It is worth to note that both solutions lead S = T25Z. This equality was also found in BSFT with 
the boundary interaction (fT9|) [2]. Since the difference between S and T25Z was given by dZ/dhifj,, this 
equality implies that the solutions correspond to the conformal fixed points of the world sheet theory. 

8 It is interesting that the second solution T = 00 lets all the other equations of motion for be 
trivially satisfied, for any value of {u^}. This is a generalization of the fact that at the tachyon vacuum 
any constant field strength of the massless gauge field is a degenerate solution. Probably this is related 
to the fact that there is no open string excitation at the tachyon vacuum. 
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To find nontrivial solutions with massive state condensation, we solve this df jdu = 
numerically by truncating the fields as 



u k = for k > k c . (38) 

The fact that the variation of f(u) with respect to u k consists of only u k and Uk±i implies 
df /diik>k c +i = 0. Hence the nontrivial equations we have to solve are 

df 

— - = , for k < k c + 1 . (39) 

OUk 

In general, there is no solution since the number of the equation is k c + 1 and is bigger 
than the number of degrees of freedom k c . We first neglect the stationary condition with 
respect to Uk c +i and find the numerical solution of 

df 

EOM(£; c ) : — — — — , for k < k c . (40) 

Let v kc,s be the solution of EOM(/c c ), where s is the natural number labeling the solutions. 
We sort solutions in ascending order in their values of /, i.e. f(v kc ' Sl ) < f(v kc,S2 ) for any 
choice of si < S2 in the set {s}. If df /du kc +i is sufficiently small at u — v kc ' s , we can 
regard v kc,s as the approximate solution of the whole system of equations (|39|) . 



3.2 Numerical solutions 

Let us show numerical solutions. We shall also explain how we can make the solution 
accurate, in spite of the introduction of the effective cut-off k c . We find a large number 
of solutions having different energies. 

Since the energy is written in terms of the function / as in (136|) . first we study the 
numerical solutions in terms of the values of /, which can help the analysis easier. 

We solve EOM(/c c ) numerically for k c — 2, 3, ■ • • ,6. The values of f(u) and \df /duk c +i\ 
for the real solutions are plotted in Fig{TJ The number of real solutions are 3, 5, 11, 21, 44 
for k c = 2, 3, 4, 5, 6 respectively. So the number of solutions diverges rapidly as we increase 

kE 

At this stage, df /duk c +\ is not so small and is about the same magnitude as f{u) 
itself. So, we cannot tell that these numerical solutions solve the full equation of motion. 
We shall improve the situation below. 

9 From this data, approximately, the number of the solutions N(k c ) increases by a factor of 2 when 
we increase k c by 1, so one can approximate N(k c ) oc 2 ka . However, the data can be fit well also with 
N(k c ) oc k^. So we do not conclude whether the growth of the number of solutions is exponential or 
power-law-like. See the discussions at fig. [7] 
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Figure 1: The left and the right graph show the values of f(v' s ) and 
k c = 2, 3, 4, 5, 6 respectively. 



I 9f 



for 



For k c > 7, the computatios become much more complicated. Hence we take an 
alternative approach. We begin with the s-th solution of EOM(/c c ), v kc,s . Using the 
Newton's method with the initial value 

u k = v k k c,s for k <k c , u kc+x = , (41) 

we obtain the solution of EOM(fc c + 1) which we call yj hc ' s,kc+1 . In the same way, we can 
find the solution of EOM(/c c + 2) by applying the Newton's method with the initial value 

Uk = w M+i for k <k c + l, u kc+2 = 0. (42) 

By iterating this procedure, we can obtain the solution of EOM(L) for any L > A^l^l 

The above procedure is shown to work well as an approximation, and provides a con- 
vergent solution solving the full equation of motion. Let w kc,s ' L be the solution of EOM(L) 
found by the iteration procedure which begins with v kc,s , the solution of EOM(fc c ). We 
found that as a function of L, f(w kc,s ' L ) converges to a nonzero finite value as L — > oo, 
while df/duL + i tends to vanish as 0(1/ L 2 ). In Fig|2l we show an example (among 
many solutions) of the plots of / and \df /du L+1 \ for k c = 6, s = 5 as a function of 1/L and 
1/L 2 respectively. In this example, / ~ 0.0071 and \df/du L+1 \ ~ 0.000017 at L = 100. 
Therefore we claim that the equation of motion is approximately satisfied. By fitting 
f(w 6,5,L ) for large L, say 30 < L < 100, to a quadratic function of 1/L, it turns out to 
approach to a non-trivial value 0.006887. 

10 It is, however, important to note that we get only a subset of all solutions of EOM(L) since EOM(L) 
is expected to have much more solutions than EOM(fc c ). 

11 It is not clear from the first principle that vanishing of df/duL+i as 0(1/ L 2 ) is enough since the 
number of fields diverges as O(L). The nontriviality of the equation of motion due to an infinite number 
of fields also appears in cubic string field theory. We leave further study of this problem to future work. 
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f(u) 



\df/du L+[ 



0.012 : 

0.011 ; 
0.010 ; 
0.009 ; 
0.008 ; 
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0.006 - 



L 1 0.000 " 



Figure 2: The plots of f(w G,5,L ) (the left graph) and | 
functions of 1/L and 1/L 2 respectively. 



w 



6,5,ZA 



(the right graph) as 



For any k c and s within 2 < k c < 6, we extrapolate the values of / at L — oo, named 
foo{k c ,s), from the data for 30 < L < 100. Here the label s of the solution is chosen in 
the same manner as before, foo(k c , si) < foo{k c , S2) for any Si < s 2 . The result is shown 
in FigH 

The horizontal lines reveal that solutions of EOM(4) include ones of EOM(2) and 
EOM(3) and are included in ones of EOM(5) and EOM(6). For this reason, we expect 
that {w kc ' s,L=00 } with fixed k c forms a subset of solutions of EOM(/c c = 00). 

3.3 Distribution of the energies of the solutions 

Based on the above results, we study the energies of the solutions. We shall concentrate 
on Lorentz-invariant solutions for simplicity in this subsection. Recall that the potential 
energy is written in terms of f(u) as (|33l) . The Lorentz invariant solutions are given by 
choosing = u for all /i. Now the energy in units of T25V26 is given by e -26 -^^ and so 
the extrapolated value at L = 00 is e - 26 /°°( fc c,s)_ plots of the energies are given in 
FigH 

We notice the following interesting features of the distribution of the energies of the 
solutions: 

(1) Almost uniform distribution of the energy spectrum. 

Although we have found a lot of solutions, the energy of those solutions do not 
overlap with each other, while tend to be uniformly distributed. This can be seen 
in FigJH as linear profiles of each dotted sector. 

(2) Energies below the D-brane tension. 
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f(u) at L=oo 

0.06" 
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kc=2 


■ 


kc=3 


♦ 


kc=4 


▲ 


kc=5 


▼ 


kc=6 



Figure 3: foo(k c ,s) extrapolated from the data for 30 < L < 100. Here we draw the 
horizontal lines in accordance with the values of /oo(4, s). 



All the solutions we found have energies which take the values between the D-brane 
tension (normalized as 1.0 in FigH|) and the tachyon vacuum (no D-brane). 

(3) Solution with a very small energy density. 

The lowest value of the energies among a fixed set decreases as k c increases. It 
suggests that bringing up k c further would lower the lowest value of the possible 
range of the energy distribution. 

(4) The presence of the "desert". 

As seen from the plots in Fig. |3l there exists a "desert" in 0.25 < foo < 0.4 where 
there appears no solution up to k c = 6. The desert, however, gets narrow as k c 
increases. For this reason, we expect that this desert is an artifact of finite k c . 

These features on energies are interesting, besides the surprising fact that we have 
obtained a large number of solutions in string field theory. The large set of the solutions 
would serve as a "string landscape." 

It requires a possible interpretation in terms of string theory. In the next section, we 
provide a possible interpretation of the solutions, with a detailed analysis of the intervals 
of the energies and the degeneracies of the solutions. 

4 An interpretation of the solutions: closed strings 

Closed strings, which should live at the tachyon vacuum as physical excitations, remain a 
mystery in string field theory. There should exist closed string excitations, in particular 
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Figure 4: The plots of the energies for Lorentz invariant solutions in units of T25V26 
extrapolated from the data for 30 < L < 100. For every k c , there is the solution whose 
energy is 1, or T 2 ^V2q- These are the trivial solutions u = 0. 



at the tachyon vacuum where the open string degrees of freedom should go away along 
the disappearance of the original D-branes a la Sen's conjecture. Now, a first look at the 
energy plot of Fig. @] would suggest a set of closed string states. The reasons are almost 
obvious: the plots appear to have a uniformly quantized energy levels (the feature (1) in 
the above list), and are below the D-brane tension, raising up from the tachyon vacuum 
(the feature (2)). 

In this section, we first present a detailed analysis on how the solutions are distributed, 
and the existence of the degeneracy of the solutions. Then we discuss that both may serve 
as indirect evidence for our interpretation of our solutions as closed string states, but the 
true connection to the closed strings states in the tachyon vacuum is yet to be unraveled. 

4.1 Properties of the solutions 
4.1.1 Uniformity of the energy intervals 

We first analyze the distribution of the function / for the solutions. To investigate the 
distribution of foo(k c , s), we evaluate the intervals of pairs of foo(k c , s): 

goo(k c , si, s 2 ) = — - — (foo(k c , s 2 ) - foo(K, si)) , (43) 

S2 ~ Si 

for any pair s± < S2- Henceforth, we don't take into account the pairs separated by 
the desert. The result is shown in Figj5j It turns out that the intervals get smaller as 
k c increases. This is the outcome of the increase of the number of the solutions. It is 
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Figure 5: The left graph shows the intervals of any paris of foo(k c , s) next to each other, 
namely, <7oo(^c> s > s + 1); while goo(k c , s±, s 2 ) for general s± < s 2 are shown in the right 
graph. 
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Figure 6: The energy intervals in units of T 25 V2 6 . The left graph shows E(k c ,s,s + 1), 
the energy intervals of any paris next to each other, while the intervals of any pairs are 
shown in the right graph. 



important to note that the intervals become almost identical to each other. This implies 
that the distribution of foo becomes uniform. 

To see the intervals of the energies we define, in the same manner, 

E(k c , si, s 2 ) = — - — ( e - 26 ^( fc -si) _ e -26 fx (k c , S2 h , 44 ) 

S2 ~ 

where we don't take into account the pairs separated by the desert. Again, the distribution 
of the energies is seemingly uniform except the desert. See Figj6) 

So, as expected from a brief look at the energy distribution of Fig. HJ it is indeed 
the case that the intervals between energy levels for different solutions are approximately 
identical to each other, and the quantized distribution of the energy levels are uniform. 

The fact that the distribution of the energies is uniform appears to be inconsistent 
with the uniformness of the distribution of f^. For the solutions we found, however, the 
values of are small so that e~ 26 ^°° ~ 1 — 26/qo and so there is not much numerical 
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Figure 7: The ratio of the averaged intervals: the averaged intervals for k c divided by 
the one for k c + 1. The left graph shows the ratio for / and the right graphs shows the 
one for the energy. 



difference between the intervals of foo and the ones of the energies. At this stage, we can't 
precisely conclude whether / or some function of / like e -26 ^ has a uniform distribution. 
To obtain more precise prediction, we need to extend our computations to higher k c where 
we might get solutions with a sufficiently large value of / to investigate which intervals 
are uniform. 

The averaged interval is smaller for larger values of k c . FigJT] shows that the averaged 
interval can be fit well as a function 2~ kc . This is consistent with the fact that the number 
of solutions grow as 2 kc as stated earlier, assuming that finally in the k c — > oo limit the 
energy levels are uniformly distributed between the D-brane tension and the tachyon 
vacuum. 

4.1.2 Degeneracy among the solutions 

Previously, we have considered only a solution with Lorentz invariance; we took all u M to 
be equal to each other. In general, however, by choosing different solutions of f(u) for 
different components of u^, we obtain solutions breaking Lorentz invariance. Apparently, 
the energy spectrum degenerates since the equation of motion is still satisfied when we 
permute the components of « M . This is an exact symmetry among solutions. 

Suppose we decompose 26 into two subsets. The I components in the first subset 
take m m = 0, while the other 26 — / components take a single solution = u%. Obviously, 
we can find 26!/(26 — solutions with degenerate energy, from how we arrange the 
subsets. So, allowing a breaking of the Lorentz symmetry provides a huge degeneracy in 
energies of the solutions. Note that we can decompose 26 components further into a large 
number of subsets. 
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There is an additional form of an approximate degeneracy which comes with our 
numerical finding mentioned earlier, the approximate uniformity of the distribution of 
the solutions in the foo space. Let us formally express solutions of EOM(/c c = oo) as 
{/oo(oo, s)}. The uniformity of the distribution implies 



where f\ = /oo(oo, 1). It turns out that there are two types of solutions which have the 
same energy e~ 2 ^: 



Therefore, the energy spectrum degenerates. This degeneracy adds up on the exact degen- 
eracy explained above. So, in total, we would have a huge degeneracy in energy, among 
our solutions. 

4.2 An interpretation: closed string states? 

Let us discuss whether the properties of our BSFT solutions may allow an interpretation 
as closed string excitations at the tachyon vacuum. As we will see, it is not conclusive. 

First, let us see the uniform distribution of the energy of the solutions. One should 
note that as the D-brane tension is inversely proportional to the string coupling constant 
g s , perturbative spectrum should appear infinitely dense if we normalize the D-brane 
tension to be the unity and take the perturbative string limit g s — >■ 0. So, the increase 
of the number of the energy levels for larger k c is consistent with the interpretation that 
those are some perturbative excitations of string theory. 

The uniformity first looks as closed string energy levels, but note that our energy is the 
bulk energy of the BSFT solutions, while in closed string excited states what is uniform 
is the closed string hamiltonian as a single-body problem. Since we do not know why 
closed strings should be homogeneously distributed in space, we lack a direct connection 
between closed strings and our BSFT solutions. 

The standard closed string states have a huge degeneracy among states sharing the 
same energy. As we saw above, we have found a similar degeneracy in our spectrum 
of the energies of the solutions. ( 1461) is reminiscent of the closed string spectrum. In 
fact, we can reproduce a part of the closed string degeneracies by our solutions, and the 
representations under Lorentz transformation of a part of closed string spectrum by the 
following simple rule of a replacement: f{u^) = nf\ — > At n = at n at n . For example, two 
solutions given in ( l46l) correspond to A^A^O) and At 2 \0), respectively. 




for s = 0, 1, ... , 



(45) 




(46) 
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There are, however, two important discrepancies: First, in our solutions, there seems to 
be no solution which amounts to the closed string state A^A^O). Second, more notably, 
the degeneracy increases as the energy decreases, in contrast to the closed string spectrum, 
since f(u) contributes to the energy in the form of e~^ u \ The latter problem would be 
complicatedly related to the infinitely large density of energy levels in the k c — > oo limit, 
which may require more exploration of the solution space at larger k c . 

In sum, although the uniformity and the degeneracy are quite suggestive, they are not 
sufficient to claim that our solutions are closed string states. It requires further study for 
a conclusive interpretation. 

5 Discussion 

In the present paper, we have solved the equations of motion derived from the BSFT action 
associated with general quadratic boundary operators. By means of numerical analysis, we 
have found a large number of solutions whose energies are uniformly distributed between 
the energy of the tachyon vacuum and the D-brane tension. As the quadratic boundary 
interactions give a free worldsheet theory, the solutions we obtained are solutions of the 
full string theory. 

As we have discussed in sec. HI our solutions are possibly related to non-trivial closed 
string excitations at the tachyon vacuum (alternatively called the closed string vacuum). 
In fact, it was discussed in [35] that the non-local open string background implements 
shifts in the closed string background. By extending our analysis to larger k c , we may 
make a progress on the interpretation of the solutions. 

Our solutions in BSFT do not have any counterparts among known solutions in CSFT. 
Although there are some suggestions on a possible relations between the two SFT's [36] . 
it is difficult to see how our solutions may be mapped to CSFT. It would be interesting if 
one can construct CSFT solutions sharing the properties with our solutions. In addition, 
to gain insight on what our solutions mean, it is important to know physical excitations 
around the found solutions. Whether the closed string excitations actually can be iden- 
tified would be the key point. Furthermore, an analogue of rolling tachyon solutions in 
BSFT [371 EH EH], and possible deformations of solutions in BSFT [ID] would gain more 
insight. One of our Lorentz- violating solution resembles tachyon matter solutions [4Tj . 

While we worked with the non-local boundary operators in this paper, the BSFT 
action associated with the "local" boundary operators can be obtained by setting u^ u {9) = 
^2 s r= o^r^ r \d) . Due to the contact divergence we need to introduce the short-distance cut- 
off e and renormalize the coupling a as a' = a + Aa(t, e) . With a proposed form of the 

12 A different form of the renormalization using the so-called ^-function renormalization is discussed in 
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counter term Aa in [20], the BSFT action S s (to, . . . ,t s ) turns out to manifestly depend 
on s: S s (t , . . . ,t s — > 0) ^ S s _i(t , . . . , Since the counter term Aa just shifts a, 

after we integrate out a, the action is independent of the choice of Aa. We found that 
this action again depends on s in the above sense. The physical interpretation of this fact 
is not clear and this is why we didn't work with the local boundary interactions. 

In the introduction, we mentioned that massive mode condensation would be related to 
multiple-D-brane solutions in SFT. Although in sec. H]we discussed our solutions may be 
interpreted as closed string states, they may still allow another interpretation as multiple- 
D-branes. We have not obtained a solution with energy larger than the original D25-brane, 
it would not mean that our solutions are not multiple-D-brane solution. The reason is 
that since we are working in bosonic string theory, multiple-D-branes are unstable and 
would form a bound state whose energy may be much smaller than just the multiple of the 
D-brane tension. This kind of question can be answered only in superstring field theory, 
as superstring should have stable multiple-D-branes with energy protected by the BPS 
property, and further study is necessary. It is, however, nontrivial to extend our analysis 
to superstring. In fact to find the action is even more involved since the associated 
boundary operators are no longer quadratic in the matter operator X M contrary to the 
bosonic string (see [13] for an explicit treatment of massive states in super BSFT). 

At a glance, our main result of figJH resembles a band structure of electrons in mate- 
rials. The fermion band structure is related to matrix models, and indeed some matrix 
models represent tachyons and unstable D-branes (see for example [HI H5]). Suppose 
our solutions with the massive field condensation are bound states of multiple unstable 
D-branes discussed above, then it is natural that matrix models appear as a low energy 
description of the multiple D-branes. The excitations of the matrix models may look like 
a band structure. In general, the number of the excitations of the matrix model increases 
as the rank of the matrix N increases. In this sense, k c plays a similar role as N. It is 
interesting to seek the counterpart of the "desert" with given k c in the matrix models 
with finite N. 

From the cosmological point of view, the distribution of infinitely many solutions is 
suggestive of the so-called landscape, which was found in superstring theory. Our solutions 
would be called as a landscape of BSFT. Among its peculiar properties, it is intriguing 
that the lowest value of the energy decreases when we increase k c . If we take a large enough 
k c , one may have a BSFT solution with very small cosmological constant. In secJH we 
also commented that we have solutions breaking Lorentz invariance. This fact suggests 

02- 

13 It is important to note again that to find the excitations around our solutions is important. There 
would be tachyonic excitations. 
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that generic non-perturbative vacua may spontaneously break the Lorentz symmetry. It 
reminds us of the original motivation of exploring non-perturbative vacua in string field 
theories in late 80's: the spontaneous Lorentz and CPT violation [2TJ [221 122]- It would 
be interesting to further explore our Lorentz-symmetry-violating (metastable) vacua and 
its relevance to the important role in the baryon asymmetry of the universe. 
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